Modeling the Unconventional Superconducting Properties of Expanded A3C60 
Fullerides 



Massimo Capone/'^ Michele Fabrizio,-''^ Claudio Castellani/ and Erio Tosatti^''* 
^SMC, CNR-INFM and Dipartimento di Fisica, Universita "La Sapienza", P.le Aldo Mom 2, 
1-00185, Roma, Italy 

^ISC-CNR, Via del Taurini 19, 1-00185 Roma, Italy 

^International School for Advanced Studies (SISSA), and CNR/INFM DEMOCRITOS 
National Simulation Center, Via Beirut 2-4, 1-34014 Trieste.ltaly 

''The Abdus Salam International Center for Theoretical Physics (ICTP), P.O.Box 586, 
1-34014 Trieste, Italy 

(Dated: September 4, 2008) 

The trivalent alkali fullerides solids of generic composition A3C6O1 where Cgo is the fuUerene 
molecule and A = K, Rb, and Cs, are a well established family of molecular superconductors. 
The superconductive electron pairing is of regular s-wave symmetry and is accounted for by 
conventional coupling of electrons to phonons, in particular by well understood Jahn Teller 
intramolecular Ceo vibrations. A source of renewed interest in these systems are alarming 
indications of strong electron-electron repulsion phenomena, which emerged especially in 
compounds where the Ceo-Ceo distance is expanded, by either a large cation size or by other 
chemical or physical means. Several examples are now known where this kind of expansion, while 
leading to a high superconducting temperature at first, gradually or suddenly causes a decline 
of superconductivity and its eventual disappearance in favor of a Mott insulating state. This 
kind of insulating state is the hallmark of strong electron correlations in cuprate and organic 
superconductors, and its appearance suggests that fullerides might also be members at large of 
that family. 

Our approach to the fullerides is theoretical, and based on the solution of a Hubbard type 
model, where electrons hop between molecular sites. We take advantage of the fact that, unlike 
models for the strongly correlated cuprates, still under debate, in a Hubbard model of fullerides 
all the important electron correlations occur within the molecular site, efficiently soluble in 
the Dynamical Mean Field Theory (DMFT) approximation. DMFT solutions confirm that 
superconductivity in this model fuUeride, although of s-wave symmetry rather than d-wave, 
shares many of the properties that are characteristic of high Tc cuprates. The calculations are 
heavy; and while our working model is several years old, the new results we present in this 
Colloquium pertain to the most interesting case of three electrons per Cgo molecule, appropriate 
to AaCeo, and have only been possible recently thanks to a stronger computational effort. 

We have calculated the zero temperature phase diagram as a function of the ratio of intra- molecular 
repulsion parameter U over the electron bandwidth W, the increase of U/W representing the main 
effect of lattice expansion. The phase diagram is close to that of actual materials, with a dome 
shaped superconducting order parameter region preceding the Mott transition for increasing cell 
volume. Unconventional properties of expanded fuUeride superconductors predicted by this model 
include: (i) an energy pseudogap in the normal phase; (ii) a gain of electron kinetic energy and 
of conducting Drude weight at the onset of superconductivity, as in high Tc cuprates; (iii) a spin 
susceptibility and a specific heat behavior that is not drastically different from a regular phonon 
superconductor, despite strong correlations; (iv) the emergence of more than one energy scale 
governing the renormalized single particle dispersion, electronic entropy and the specific heat 
jump. These predictions, which if confirmed should establish fullerides, especially the expanded 
ones, as members of the wider family of strongly correlated superconductors, are discussed in the 
light of existing and foreseeable experiments. 
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this theory, electron-electron repulsion merely renormal- 
izes the electron-phonon parameters, lowering the critical 
temperature Tc rather than enhancing it. 

The trivalent alkali ful lerides superconductors, e xcel- 



Dis covered by K amerlingh Onnes nearly a century 
ago ( Onnesl . Il91lh . and first explained microscopi- 
cally back in 1957 in terms of elect ron pairing by 
Barde en, Cooper and Schrieffer (BCS) (jBardeen et all . 
1957al lbl). superconductivity is still surprisingly a' la page. 
On one hand, superconductivity is being constantly dis- 
covered in an ever increasing variety of solid state com- 
pounds. On the other hand, it appears more and more 
difficult to use basically the same st andard theo ry, es- 
sentially BCS and its extensions (e.g. iParksl Il969l chap- 
ters 10 and 11) to account for all of them. In this stan- 
dard, conventional theory, superconductivity arises from 
the condensation of electron pairs, the two electrons usu- 
ally bound in a pair state of s-wave symmetry and held 
together by exchange of lattice phonons. The Coulomb 
repulsion between the two electrons opposes pair forma- 
tion, but it does not suppress superconductivity because 
screening makes it weak enough. 

The surprisingly favorable effect of repulsive electron 
correlations on superconductivity found in some sys- 
tems, particularly in high- Tc superconducting cuprates, 
where electron-electron repulsion is dominant, remains 
a standing puzzle. An immense amount of experimen- 
tal and theoretical work has accumulated over the last 
two decades in the attempt t o understand these phenom- 
ena, s ee e.g. the review by iBennemann and KettersonI 
( 20081 ). Actually, cuprates are but the most spectacular 
members of a wider class of strongly correlated supercon- 
ductors incl uding heavy fermion and organic molecular 
compounds ( Bennemann and KettersonI |2008[ ). systems 
for which there is no really comprehensive theory either. 
Among other factors, theoretical efforts have been ham- 
pered by the general inter-site nature of electron inter- 
actions and correlations in many of these systems, a fact 
that poses large technical difficulties. In this light, iden- 
tifying a superconductor family where correlations are at 
the same time strong, simple, and on-site is welcome. 

A more crucial element is one of perspective. It has 
been a widespread prejudice to distinguish between su- 
perconductors where (as in BCS theory) pairing of elec- 
trons takes place in the s-wave channel and is mediated 
by phonons, from those where the mechanism may be 
electronic and not phononic, and where pairing instead 
takes place in the d-wave channel. Whereas it is widely 
held that strong repulsive corr elations are essen tial to su- 
perconductivity in the latter ( Andersonl . Il987t ). they are 
not considered crucial in the former. The conventional 
BCS scenario an d its extensions, namely the Migdal- 
Elias hberg theory (|Enashberel[l960l:lMigda]| . ll958HPark£l 
19691 ) - a controlled approximation valid when the typ- 
ical phonon frequency is much smaller than the Fermi- 
energy - are more or less automatically accepted, and 
used to account for the superconducting properties. In 



lently re v iewed , e.g., by [Cunnarssor 
IRamired (I1994D . are among the systems where this con- 
ventional logic seemingly applies. Fullerides are solid 
state compounds of generic composition A3C60 , where 
Cgo is the fuUerene molecule ( Gunnarssonl . |20M ) and A 
= K, Rb, and Cs are alkali cations. The three alkalis 
donate a total of n=3 electrons to each fuUerene, half fill- 
ing its threefold degenerate tiu molecular level. Electron 
hopping between first neighboring fuUerenes gives rise to 
a metal, where conduction is restricted to the three nar- 
row ii„-derived bands, w ith a t o tal energy bandwidth 
of no more than 0.6 eV (jErwinl . Il993t ISatpathv et~d\ . 
I1992D . Metallic fullerides are generally superconducting, 
with critical temperatures Tc reaching ~ 40 K, depend- 
ing on various factors. An empirically important factor 
appears to be the cell volume. When the fuUeride lattice 
is chemically expanded, by either increasing cation size 
or by insertion of neutral molecules, or else physically 
expanded by removing pressure, Tc undergoes a definite 
and systematic change. It rises initia lly with a good cor- 
relation with the Ca o-Cfio distance ( Gunnarssonl Il997t 
lYildirim oil . I1995D . Further expansion however causes 
Tc to drop, ending eventually, through a first order tran- 
sition, in an insulating state, as we shall discuss later. 

A wealth of evidence indicates that superconducting 
pairing in fullerides is phononic and that the relevant 
phonons are the stiff intra-molecular Hg vibrations of the 
Cgo molec ule, Jahn- Tel l er cou pled to the tiu conduction 
electrons ( Gunnarssonl [l993)- Further support to the 
apparent BCS nature of superconductivity in fullerides 
comes from specific-heat jumps that scale linearly with Tg 
in ag r eement with BCS theorv dBurkhart and Meingast , 
119961 : iKortan et all . Il992t [Ramired, llSM), as well as a 
regular (i.e., not exceeding ly high ) normal phase mag- 
netic susceptibility. (Kortan et all Il992t [Robert et all . 
Il998[ ) S uperconduct i ng ene rgy gaps are less clearly de- 
fined, ( Gunnarssonl Il997[ ) the ga p ratio 2A /Tc ^ 

3.4 - 4.2 in K sCrq and RbaCeo (|Gunnarssoil Il997t 
iRamirej . 1 19941 ). but not far from the BCS value of 3.53. 
These elements suggest viewing the A3 Ceo c ompounds 
as we akly correlated Fermi-liquid conductors ( Ramire3 . 
I1994D . though with unusually narrow electron bands, 
with a large effective mass roughly three free electron 
masses ( Robert et all . \l99§t) . Even the observed decrease 
of Tc under applied pressure in KaCgo and RbaCeo is in 
qualitative agreement with an increasing bandwidth and 
decreasing density of states at the Fermi level, which fur- 
ther supports a standard BCS picture. 

These reassuring, conventional looking facts are 
however contrasted by a number of conflicting ele- 
ments that are strong enough to cast serious doubts 
on the general applicability of the BCS scenario 
to superconductors in this family. These elements 
are especially apparent in the more expanded ful- 
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lerides inclu ding: (NH^lxNaK^Cfin (|Ricc6 et all l2003l ) 
( Durand et ai] . 1200% 



and in the alloys 



alloy s 
1200(1 ). 



and LigCeo 

Cs3_^Ka;C6o and Cs3_a;Rba;C6o (jPahlke et al 
For these expanded compounds Tc decreases upon ex- 
pansion, contrary to BCS theory. The electron density of 
states extracted by NMR Knight shift is at the same time 
an increasing function of lattice parameter, smoothly 
connectin g with that o f the u nalloyed compounds under 
pressure (jPahlke et aZ.I . l2000l ). Within BCS theory, that 
increase should lead to a rise of and not to a drop 
as observed. The same unconventional behavior is ob- 
served in CsaCgo, the fuUcride compound with the high- 
est T g ~ 40 K attained under pressure ( Palstra et al\ . 
I1995D . A novel A15 superconducting phase of CsaCeo 
with expand ed structure ha s very recently been syn- 
thesized (|Ganin et al\ . l2008l ) corresponding to a body- 
centered cubic arrangement of fullerenes. Superconduc- 
tivity emerges under pressure through a first order non- 
structural transition at 4 Kbar. The critical temperature 
Tc first increases with pressure, reaching a dome-shaped 
maximum of 38 K around 7 Kbar, above which drops. 
Since no structural changes are observed under pressure, 
the appearance of superconductivity as well as the dome- 
shaped Tc vs. pressure be havior must be ascribed solely 
to the volume contraction ( Ganin et Hl,[2008). This non- 
monotonic behavior of Tc with pressure finds no apparent 
explanation within the conventional theory. 

The basic and striking anomaly of expanded fuUerides 
occurs in the compounds with the largest inter-molecule 
distances. In these materials a relatively modest ad- 
ditional lattice expansion (and minor change of sym- 
metry due to intercalated ammonia) is enough to dra- 
matically turn them from metallic and superconduct- 
ing to antiferromagnetic and ins ulating ( Durand et al\ . 
120031 : llwasa and Takenobul . l2003l ).^ With temperature, 
antiferromagnetic order in the ammoniated compound 
NH3K3C60 changes to paramagnetic di sorder at a Nee l 
temp erature slightly above ~ 40 K ( Prassides et al\ . 
I1999I ). Even above the Neel temperature, the microwave 
conductivity in NH3K3C60 remains several orders of mag- 
nitude below that of K3C60 (|Kitano et al\ . l2002t ). testi- 
fying the Mott insulator (correlation driven) nature of 
the insulating phase. Electrons in a lattice give rise to 
a Mott insulating state when electron electron repulsion 
stops their free propagation and the lattice appears as a 
collection of molecular ions. Correlations lead to an en- 
ergy gap in their spectrum, even if their number density 
per c ell is odd, instead of even as in regular band insula- 
tors ( Mottl . I199OI ). Proximity of a Mott insulator phase 
in fuU erides had long been advoca t ed in different con- 



texts (iBaskaran and Tosatti . 1991: Chakra vartv et 



ty 

I1991I: IChakravartv and Kivelson. . ,1991; .Lof et all Il99a ) 




^ We expect that below the super conducting pressure of 4 kbar, 
the new compound A15 CsaCeo jGanin ef am200§) . yet to be 
characterized in this respect, should also be an insulating anti- 
ferromagnet. 



FIG. 1 Schematic representation of a planar projection of 
the crystal structure of NHsKsCeo- Green dots are potas- 
sium atoms, blue dots surrounded by three grey dots are NH3 
molecules, and Cgo molecules are shown according to their 
actual spatial orientation. (Courtesy of Kosmas Prassides). 



but was not taken seriously by the community prior to 
this data. Superconductivity next to a Mott insulating 
phase as a function of doping or volume change is the 
hallmark of strong correlations in high temperature su- 
perconducting cuprates and organics. One is thus nat- 
urally led to inquire whether superconductivity in ex- 
panded fuUerides might, despite the obvious differences, 
and despite the phononic mechanism, be somehow re- 
lated to strong correlations. Our contention is that it is 
indeed closely related, as outlined in the following. 



II. ELECTRON INTERACTIONS IN FULLERIDES 

The electrons donated by the alkalis to the Cgo 
molecule enter the threefold degenerate tiu former low- 
est unoccupied molecular orbital (LUMO). In a degen- 
erate molecular orbital, electrons interact through a va- 
riety of mechanisms. The first is overall Coulomb re- 
pulsion, which we will associate later below with the 
Hubbard parameter U . The second is Coulomb ex- 
change energy, minimized when the molecular state 
has the highest total spin, and the highest total or- 
bital a ngular momentum c o mpat ible with it (Hund's 
rules) (|Landau and Lifshit j . Il958l ). The third is the 
Jahn- Teller (JT) interaction, caused by coupling of the 
electr on levels to symmetr y lowe ring molecular distor- 
tions ( Landau and Lifshitd . Il958l chap. 13). Contrary 
to Hund's rules, in a JT distorted molecule the ground 
state maximizes double occupancy of levels, thus favor- 
ing low total spin instead of high spin in the isolated 
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molecular ion. In molecular Cgo", the strength of these 
interactions has been evaluated in the past, and the JT 
strength has been esti mated to prevail n arrowly over 
Hund's rule exchange i Liiders et all |2002( ). This nar- 
row balance favors a low-spin ground state, with a rel- 
atively small "spin gap"- the energy between the low 
spin ground state and the lowest high-spin excit ed state , 
expected to be of the or der of 0.1 eV ( Capone et all . 
I2OOII : iLtiders et al. In agreement with this ex- 

pectation, local moments indicate that in antiferromag- 
netic Mott insulating NH.'^ K.^Cfin the Cfin(3— ) sit es are 
in a low-spin state, S =1/2 (jPrassides et all . [l99^ . their 
high spin state 5 = 3/2 lying about 100 meV higher 
in energy. A low-spin qualifies the overexpanded ful- 
lerides as "Mott- Jahn- Teller" insulators - that is, Mott 
insulators whose sites are in a JT stabilized low-spin 
state ( as opposed to a Hun d 's rul e stabilized high-spin 
state) (|Fabrizio and Tosattil. 119971 ). Under hydrostatic 
pressure, NH3K3C60 undergoes a transition to a metallic 
state wh ere superconduc t ivity re-emerges with a rather 
large (jPrassides et a/.l . ll999D . It is important to note 
that this superconducting phase still belongs to the "ex- 
panded" family, as signaled by the fact the Tc here in- 
creases furthe r with increasing pressure (reaching 28 K 
at 14.8 kbar (|Zhou et al\ . [1991)) at variance with non- 
expanded fullerides, where Tc drops under pressure. 

A lateral, but relevant additional element comes from 
tetravalent compounds A4C60, which are insulators or 
near insulators. By comparison with the trivalent ful- 
lerides, the slight reduction of band-energy gain per par- 
ticle caused by adding one more electron per molecule 
and by slightly changing the crystal structure is suffi- 
cient to turn the trivalent metals into tetravalent insu- 
lators even in non-expanded materials. Careful density- 
functional electronic structure calculations indicated that 
it is not possible to describe the tetravalent compounds 
such as K4C6n as statically di storted Jahn- Teller band in- 
sulators ( Capone et a/.l . [2000h . A static JT distortion and 
the associa ted orthorhombic state is actu ally present only 
in CS4C60 (jPahlke and Rosseins ky, 2002 V and in R b4C6o 
above a critical pressure (iHua and Stephena. 2006 ), while 



it never shows up in K4C60 (|Hua and Stephensl. 120061) 
(with the exception of monolayers, see IWachowiak et al\ 
( 20051 )). The persistence of insulating or near insulat- 
ing behavior and the recovery of molecular symmetry 
observed in the high temperature phase of tetravalent 
fullerides suggests tha t these compounds too are Mott- 
Jahn - Teller insulators ((Capone et al. l. l2000l : lKlupp et all 
I2OO6I : iKnupfer and Finkl . I1997D 7 like the overexpanded 



trivalent materials. The dynamic JT effect in each Cg,^ 
ion associated with Mott localization of carriers is cru- 
cial in explaining the low spin ground state and the spin 
gap of A4C60, exactly as for the expanded trivalent com- 
pounds. 

From the above discussion one might be tempted to 
conclude that strong correlations play a role only in 
tetravalent and expanded trivalent compounds, while 
face-centered cubic (f.c.c.) K3C60 and RbaCgo, where su- 



perconductivity was originally discovered, could still be 
viewed as weakly correlated systems, and as BCS type 
superconductors. We do not believe in this conclusion. 
A final, independent and strongly unconventional signal 
is provided by NMR. In fact, NMR data show direct ev- 
idence of a spin gap of order 0.1 eV, appearing as an 
anomalous activated increase of inverse relaxation time. 
Most likely this gap between a low spin ground state 
and a high spin excited state reflects the multiplet be- 
havior of the localized Cg,^ molecular ion. It shows up 
ubiquitously in all alkali doped f ullerides, including s u- 
perconducti i ig f.c. c. compounds ( Brouet et all [2002al lbl: 
iThier et~al\ . Il995l ). The existence of the spin gap sig- 
nifies that the magnetic response of fullerides is very 
far from Fermi-liquid behavior, which has no such fea- 
ture. Magnetically, the fullerides behave as if localized 
molecular multiplet excitations coexisted with delocal- 
ized propagating quasiparticles. As discussed, the recov- 
ery of molecular physics is characteristic of Mott insu- 
lators, suggesting that the fingerprint of Mott physics is 
strongly present already in the non-expanded supercon- 
ducting f.c.c. compounds. This suggests that the f.c.c. 
compounds are somehow the analog of the overdoped 
cuprates, whereas the expanded trivalent materials are 
analogous to the underdoped cuprates. The conclusion 
is that both are crucially, even if differently, influenced 
by electron correlations. 

We believe that the above elements are strong enough 
to call for a new physical picture for the whole fam- 
ily of A3 Ceo superconductors. Proximity of the Mott 
insulator strongly suggests that the anomalies of ex- 
panded fuUerene superconductors most likely originate 
from strong repulsive electron correlation in the narrow 
tiu bands. The prevalence in the Mott localized state 
of molecular physics, with its orbital degeneracy, JT ef- 
fect and intra-molecular exchange must be taken into ac- 
count along with itinerant electron band physics. Upon 
expanding the cell volume, the intermolecular hopping 
of electron weakens, whereas all the on-site correlation 
terms - Coulomb and exchange electron-electron inter- 
actions as well as molecular JT effect - are likely to be- 
come increasingly relevant. We are led to a picture where 
the Mott insulator physics of weakly coupled molecular 
ions progressively prevails over band physics for increas- 
ing lattice expansion. In particular, superconductors that 
operate in this regime are bound to deviate from the 
standard Migdal-Eliashberg scenario, the more so as the 
lattice spacing increases. To investigate that regime, we 
need to start with a broader theoretical scheme for triva- 
lent fullerides, capable of describing their behavior under 
lattice expansion and near the Mott transition. While 
that has been the scope of our work for several years, 
previous work was for pract ical technical reasons lim- 
ited to tetravalent systems ( Capone et all [200^ l200d 
I2OOII ). The study of A3C60 systems, computationally 
much heavier due to the simultaneous relevance of mag- 
netic and orbital ordering, has only now been completed, 
and we offer here an outline of the main results. 
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III. MODEL AND INTERACTIONS 

Our theoretical model of trivalent fuUerides assumes 
a lattice of molecular sites, each representing a Ceo 
molecule. The Ceo tiu threefold degenerate LUMO can 
for all purposes be treated as an atomic p level. An aver- 
age of three electrons per molecule are donated by alkali 
atoms and partially fill these orbitals, which can host up 
to six electrons. The electrons hop from site to site giv- 
ing rise to half-filled bands of width W ^ 0.6eV. On 
each site the electrons experience a Hubbard repulsion 
U ^ leV (corresponding to the Slater integral Fq oc U), 



U 2 
'Hu = ^ {n-3) , 



(1) 



together with a weaker, yet crucial, inter-orbital Hund's 
rule exchange coupling term proportional to the 
Slater integral F2. Under the sole assumption of full ro- 
tation al orbital symmetry , this exchange term takes the 
form (jCapone et al\ . l200l l200lh 



= J (^2S • S + • Lj + ^ J (n - 3)' , (2) 

where J = —Jh < 0, while n, S and L are the den- 
sity, spin and orbital angular momentum operators, de- 
fined as for p-orbitals, on the given site. This exchange 
term favors high S and L molecular multiplets, and 
has been overlooked or neglected in most treatments 
so far. The next interaction is the JT intra- molecular 
coupling of electrons in the tiu orbital to Hg intra- 
molecular vibrations. This interaction has on the con- 



trary been ver y widely discussed (s ee e.g. lAuerbach et al. 
([1994); Gunn arsson et al\ ([19951 ): iLannoo et al\ (|l99lh : 
rVarma et al\ ( 1991 ) and references therein), and we will 
not dwell too much on its details here. It acts to split 
the tiu orbital degeneracy and thus favors low spin states, 
effectively playing the opposite role to intra-molecular ex- 
change. A proper treatment of the JT coupling involves 
the dynamics of carbon nuclei in Cg^j" ions, including re- 
tardation, and has been developed bvlGun narsson et al\ 
([1995) However for expanded fuUerides close enough to 
the Mott transition, retardation is not essential and can 
be omitted. In fact, in this regime the in ter-site motion 
of qu asiparticles is severely slowed down ( Georges et all 
Il996l ) and the coherent bandwidth will eventually drop 
from W to ZW {Z < 1). When ZW falls enough to 
approach the relatively high Hp vibration frequencies o f 
fuUerene tiuj ~ 90 meV (ICapone et aLl . l200ll200ai200lh . 
quasiparticles begin to move on a comparable time scale 
with the vibrating carbo n atoms, and non-adiab atic ef- 
fects become important ( Cappelluti et all I2OOOI ) . Even 
closer to the Mott transition, the phonon dynamics be- 
comes eventually faster than inter-molecular quasiparti- 
cle hopping. In this anti-adiabatic limit, molecular vi- 
brations can be integrated away, and the resulting un- 
retarded effective JT interaction recovers identically the 
same form as Hund's rule exchange ^ except for the sign 



of J, namely Jjt > ( Capone et aLl . [200l[ ). In this hmit 
the Hund and JT intra-molecular inter-orbital interaction 
terms can be directly combined in the for m ^ with J = 
—Jh + Jjt- In Ceo, Jh ^ 0-03 — 0.1 eV ( Liiders et all 
l2002t iMartin and Ritchid. Il993l). whereas Jtt ~ 0.06 



0.12 eV (Auerbach et al . 19941 Gunnarsson ami99£ : 



'Lann oo et al . 199ll : L^ders et a7l ~ 2002t Varma et 
1991). The total result - and the only one compatible 
with s-wave superconductivity, with a spin 1/2 Mott in- 
sulator, and with a moderate spin gap near 0.1 eV - is a 
relatively weak unretarded attractio n which has the form 
of an inverted Hund's rule coupling ( Capone et aZ.I . [20021 : 
iGranath and Ostlundl . [200I . This is the approximation 
we shall adopt, keeping in mind that it is strictly valid 
only when ZW < huj, that is, close enough to the Mott 
transition. 

The lattice expansion characterizing the expanded ful- 
lerides is believed to have little effect on either Jh or 
Jjt- On the other hand, expansion will surely de- 
crease W and increase U, so it can be modeled as a 
gradual increase of U/W, reflecting both the band nar- 
rowing due to smaller overlap between molecular wave 
functions, and a decreased screening strength. As dis- 
cussed above, this Hamiltonian model, even if not really 
simple, has many body interactions that are strictly on 
site, an id eal situation fo r Dyna mical Mean-Field Theory 
(DMFT) ([Georges et a/.l . Il996[) . one of the most popular 
and powerful tools in the field of strongly correlated elec- 
tron systems that we briefiy describe in the next subsec- 
tion. 



] A. Dynamical Mean-Field Theory 



DMFT is a quantum version of classical mean-field the- 
ory, which provides an exact description of the local dy- 
namics, at the price of freezing away all spatial fiuctua- 
tions. The mean-field scheme is formulated by a mapping 
of the lattice model onto an Anderson impurity model 
(AIM) embedded in a free-elect ron Fermi "bath" sub ject 
to a self-consistency condition ( Georges et aL[ . [l99^ . In 
our model, the effective AIM is three-fold orbitally de- 
generate, with p-like levels representing the tiu orbitals, 
each hybridized with a bath. The Hamiltonian is 



n = n. 



u- 



+ Vl4a(cL,Pa, + i^.C.), 



(3) 

where T-iu and Tij are ([1]) and ^ for fermions on the 
impurity orbitals, eua are the bath energy levels labeled 
by the index k and by an orbital index a = x,y, z, Vka are 
the hybridization parameters between the bath fermions, 
created by cj.^^., and the impurity fermions, created by 
p\^. The mean- field scheme implies a self-consistency 
condition that depends on the impurity Green's func- 
tion and on the bare density of states of the original lat- 
tice. Throughout our calculations, we use an infinite- 
coordination Bethe lattice, whose density of states is 
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semicircular. This is a reasonable description of a realis- 
tic three-dimensional density of states devoid of acciden- 
tal features such as van-Hove singularities. It is moreover 
particularly convenient since it leads to a very simple and 
transparent form of the self-consistency condition. The 
Bethe lattice is the z ^ oo limit of a Cayley tree of coor- 
dination z, scaling the nearest- neighbor hopping in each 
one of the z directions as t/^/z. The resulting semicircu- 
lar density of states has bandwidth = 4t. 

The self-consistency condition requires the so-called 
Weiss field (i.e., the non- interacting Green's function of 
the AIM) 



E 

k 



lUlri — £ka 



(4) 



to be related to the local interacting Green's function, 
GaiitUn), not only through the Dyson equation for the 
AIM 

which requires the full solution of the impurity model, 
but also by the additional selfconsistency equation 



(5) 



DMFT for a given model thus amounts to solve it- 
eratively the AIM until the impurity Green's function 
satisfies Eq. ([5]). In this work, we solved the threefold 
degenerate AIM by exact diagonalization. That requires 
truncating the sum over k in Eqs. ^ and ([S]) to a fi- 
nite and relatively small number of baths N},^ so that the 
Hamiltonian can be diagonalized in the finite resulting 
Hilbert space. We generally used A^f, = 4 for each or- 
bital. Most of the results we shall present are at zero 
temperature, where we can use the Lanczos algorithm to 
calculate the Green's function without fully diagonalizing 
the Hamiltonian. The finite temperature results for the 
specific heat and its jump at Tc reported in Sec. |V] are 
the exception. They are obtained by means of the finite - 
temperature extension of Lanczos ( Capone et aZ.I . [20071 ) . 
where the thermal Green's function is expressed as a sum 
over the low-lying eigenvectors \n) and eigenvalues of 
the impurity model 



Gn 



{iuJn) 



(6) 



with 



\{n\paa\m)Y 



En 



En 



+E 



En Ei-fi 



(7) 

The Boltzmann factor in Eq. ([6]) guarantees that the 
lower the temperature, the smaller the number of ex- 
cited states that actually contribute to the Green's func- 
tion. Therefore, for sufficiently low temperature, we can 



still use the Lanczos algorithm to find the lowest energy 
eigenstates. In practice, only a limited number of states, 
20-^25, can be calculated in a reasonable time. Unfortu- 
nately, this number is insufficient to make the truncation 
error in Eq. ([6]) negligible for the present model in the 
relevant correlated regime. We estimate the systematic 
error on the Green's function to be of order of a few per 
cent, a level of accuracy that does not allow to fully deter- 
mine thermodynamic properties such as the critical tem- 
perature. Luckily though, this error affects much more 
the absolute specific heat value than its relative changes, 
including the superconducting jump in units of Tc which 
we shall discuss further below (see Sec. E]). We under- 
line that this limitation does not affect by any means the 
T — calculations. 

The above DMFT equations refer to a paramagnetic 
phase where no symmetry breaking occurs. However, in 
this work we shall be crucially interested in s-wave su- 
perconducting and in antiferromagnetic solutions, where 
symmetry is broken. Superconductivity is conveniently 
studied in the Nambu-Gor'kov representation by intro- 
ducing spinors 

^k. = ( >T \ , 

and defining accordingly the Green's function 

Gka(T) = -(TV (V'ka('^) V'L(O))) in imaginary time 
as a 2 X 2 matrix that satisfies the Dyson equation in 
Matsubara frequencies 



Gka{^to) = Gi^iiLo) + Gl^iiLu) E(zw) GU^lu), (8) 

with G^^iiuj) the non-interacting value. The single- 
particle self-energy S(ia;) is also a 2 x 2 matrix whose 
off-diagonal element A(iw), when finite, signals a super- 
conducting phase. The DMFT self-consistency can be 
written now as 

GQ^iiuJn) = il^n To+ flTs-t'^ T3 G{iuJn) T3 , (9) 

where tq and are Pauli matrices. 

Analogously, antiferromagnetism in a bipartite lattice 
is conveniently described using the spinor 



■0ka(T = 



Pkacr 
Pk+Qat 



with k in the magnetic Brillouin zone and Q the mod- 
ulation vector. This leads once more to a 2 x 2 Green's 
function and self-energy matrices related by the Dyson 
equation ([5]). Here too, a finite off-diagonal element sig- 
nals an antiferromagnetically ordered phase. The self- 
consistency equation for antiferromagnetism exploits the 
bipartite property of the lattice. Indicating one sub- 
lattice with A and the other with B, the general self- 
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consistency condition is 



(10) 



which becomes Eq. ([5]) if GaA = G^b, i-C- if the system 
is nonmagnetic. When the system is antiferromagnetic, 
then Ga-A = G-aB- Thus we can ehminate the sublattice 
B from Eq. (fTO|) . and obtain the following result for the 
self-consistency equation: 



C/q ^{iuJn)c,A = {it^n + M) - t^G{iuJn)- 



-oA- 



(11) 



This equation is valid for a Bethe lattice with nearest 
neighbor hopping, a case with perfect nesting that is 
rather exceptional in realistic antiferromagnets. A way 
to simulate imperfect nesting typical of more realistic sit- 
uations, while still taking advantage of the Bethe-lattice 
simplifications, is to add a next-nearest-neighbor hopping 
t' I z in the Cayley tree (merely a device to eliminate nest- 
ing, not meant to suggest next nearest neighbor hopping, 
small in fuUerides). In the limit z ^ oo of the Bethe lat- 
tice, the self-consistency equation becomes 



crA 



= (iw„ + /i) - G{iuJn)-aA - t' GiiUJnlii) 

For both broken-symmetry phases, if the diagonal el- 
ements of the self-energy matrix at low Matsubara fre- 
quencies follow the conventional Fermi-liquid behavior - 
which we always find to be the case. 



iagonal 



1-1 

Z 



(13) 



the actual value of the spectral gap in the single particle 
spectrum is given by Z A(0), the zero frequency anoma- 
lous (superconducting or antiferromagnetic) self-energy 
multiplied by the so-called quasiparticle weight Z. 

We end by noting that, within DMFT, one can search 
for solutions with different symmetries by simply al- 
lowing/preventing symmetry breaking order parameters, 
even in regions where the chosen phase is not the most 
stable. When two or more solutions coexist, the stable 
one is determined by an explicit energy calculation. As 
we shall discuss, for a wide range oi U/W values wc do 
find coexisting superconducting and antiferromagnetic 
solutions, the former prevailing at smaller U, the latter 
at larger U. The physical phase diagram thus exhibits a 
first-order transition between these two symmetry broken 
phases - a nonmagnetic s-wave superconductor, and an 
insulating spin 1/2 antiferromagnet - taking place when 
the respective energy curves intersect. 



B. T = Phase Diagram 

Modeling lattice expansion of fuUerides as a gradual in- 
crease of U /W, we can proceed to analyze the theoretical 
zero temperature "phase diagram" of our model obtained 



by DMFT as a function of U/W. Starting with the un- 
correlated system with U — the model initially exhibits 
straight BCS superconductivity driven by JT phonons 
(i.e., the attractive J we introduced above) with an s- 
wave {S = L = 0) order parameter^ 



(14) 



a—x,y,z 



where TV is the number of molecules and pj^g. creates 
an electron on molecule i, with spin a and in orbital 
a = x,y,z. The Fermi-liquid scattering amplitude in 
the Cooper channel, measuring the strength of the ef- 
fective attraction, is A = — lOJ/3. We note that ow- 
ing to fairly strong JT interactions ( Auerbach et all . 



Il994 iGunnarsson et all . Il995( ). if Hund's exchange Jh 
were neglected, the dimensionless JT coupling constant 
of fuUerene controlling superconductivity would be nu- 
merically very large, A = Po|^| — 0.6—1.0, where 
po (~ 2.4 eV~^) is the bare density of states per spin 
and orbital. Turning on a weak Coulomb repulsion U 
on top of that will reduce the pairing attraction in this 
regime. Perturbatively one obtains for small U that 
A — —10 J/3 + U. Since J is insensitive to expansion 
while U/W increases, this implies that in this picture, 
where Hund's rule exchange is neglected, Tc should al- 
ways decrease upon expansion, a prediction which is at 
odds with experiments. 

In fact, as anticipated, Hund's rule exchange is not 
negligible, and its effect is to introduce a substantial 
cancellation in J leading to a largely reduced effective 
coupling Xeff ~ ^ {JjT — Jh) P n- Should we interpr et 
the ubiquitous spi n gap of 0.1 eV (jBrouet et aLl . l200"2al |bl: 
iThier e^"all . [imI) as the Cg^j" molecular excitation en- 
ergy between its low spin ground state with 5 = 1/2 and 
S — for n — 3 and n — A, respectively, and its high 
spin excited state, with 5 = 3/2 and 5=1 for rt = 3 
and n = 4, respectively, we would conclude that, to be 
con sistent with our mod el where both gaps are equal to 
5 J (jCapone et all\20Q4 . the total effective inverted ex- 
change J, comprehensive of both JT and Hund's rule 
exchange, should be approximately J = 0.02 eV. In real- 
ity, the qualitative scenario we will describe is relatively 
independent of the precise value of J provided it is in- 
verted, i.e., negative and small. 

Once exchange is included, the resulting Ae// — 0.16 
is now much smaller - in fact way too small to explain 
within conventional BCS or Migdal-Eliashberg theory 
any of the observed values of Tc in fuUerides (let alone the 
non-monotonic behavio r of Tc versus the density of states 
for ex panded fuUerides ( Dahlke et a/.l . l20dol : lGanin et all 
l2008f )l. Things get worse when we increase the on-site 
Coulomb repulsion U closer to realistic values, U ^ W 



^ Since gauge symmetry is broken, we are allowed to assume Psc 
real 
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FIG. 2 Superconducting solution (triangles) and antiferro- 
magnetic solution (squares) zero-frequency anomalous self- 
energies A(0) as function of U/W. Solid symbols are used 
when the corresponding symmetry broken phase is stable, 
while open symbols when it is metastable, i.e. has lower 
energy then the other phase. The first-order transition be- 
tween the two phases is indicated by a vertical line separating 
the superconductor from the antiferromagnetic Mott insula- 
tor. Panel (b) corresponds to the case in the presence of a 
frustrating next-nearest neighbor hopping t' — 0.3t, absent in 
(a). 



and beyond. In the conventional weakly correlated pic- 
ture U would provide in the electron pairing problem a 
repulsive "Coulo mb pseudopo tential" whose bare value 
is //* = C/po - 3 (|Parksl . ll969f) . Simply comparing these 
bare values of A and /Lt* we should conclude that s-wave 
BCS superconductivity in fuUerenes is simply impossible 
(with the obvious proviso that for small U an unretarded 
treatment of JT phonon interactions is not really justi- 
fied). 

The full DMFT solution of the model for J = 0.05 W 
and increasing U/W yields the phase diagram in Fig. [21 
While confirming the above expectations for moderate 
U, it has a surprise in reserve at larger U values, where 
the Mott transition is approached. 

Fig. [3] shows the zero frequency anomalous single- 
particle self-energies calculated within DMFT for the su- 
perconducting and for the antiferromagnetic solution. At 
U = 0, the model is an s-wave BCS superconductor, 
with an exponentially small superconducting A(0). It 
is too small to be visible in the figure, since the effec- 
tive exchange- reduced A ^ 0.2 is as was said very weak. 
Beginning from zero, the increase of U first rapidly de- 
stroys the weak BCS superconductivity. The supercon- 
ducting A(0) vanishes at roughly the mean-field value 
U = 10/3 J, and above this value of U the ground state 
becomes a normal metal as expected. Upon further in- 
creasing U/W the model remains a normal metal - no 
superconductivity, no antiferromagnetism. However, the 
importance of electron correlations increases with U, as 
signaled for example in the DMFT spectral function (not 
shown) by the gradual formation of incoherent Hubbard 
bands on both sides of the Fermi level. The metallic 
character persists until a hypothetically continuous Mott 



transition eventually reached near U/W ^ 1.5, where 
Z = and the metallic character is extinguished. 

Before this happens however, s-wave superconductiv- 
ity re-enters from the normal metal state. The anomalous 
self energy A(0), proportional to the superconducting 
order parameter has, as a function of U, a bell-shaped 
behavior - a "superconducting dome" as it is called in 
cuprates - hitting a large maximum before dropping 
again. The re-entrant superconductive behavior is a clear 
realization of phonon-induced strongly correlated super- 
conductivity (SCS) (Capone et ai, 2002). The sharply 
rising order parameter edge with increasing U/W can 
in our view explain the strong rise of Tc upon lattice ex- 
pansion in non-expanded compounds, previously (and we 
believe incorrectly) attributed to a BCS-like increase of 
density of states upon band narrowing. Past the dome 
maximum, and upon increasing expansion, the SCS su- 
perconducting order parameter declines, and would even- 
tually drop to zero at the continuous metal insulator 
transition near U/W ^ 1.5. This continuous decline of 
superconductivity is preempted by a first order transi- 
tion to a lower energy antiferromagnetic Mott insulating 
phase, with order parameter 



Pa 



where 



FM 



nil), 



(15) 



is the full occupation 



-x.y,z Pi acrPi acr 

number with spin a at molecule i. The exact location of 
the superconductor-insulator transition depends on de- 
tails. For strong nesting (i'=0) it takes place even be- 
fore the superconducting dome maximum. In Fig. [3] we 
show A(0) of Fig.[2]in comparison with the spectral gaps 
ZA(0), as well as the order parameters Psc and Pafm- 
Notice that Z for the superconductor is smaller and van- 
ishes at the continuous metal insulator transition, while 
Z for the antiferromagnet is of order 1. In both cases 
the dimensionless order parameter essentially follows the 
behavior of the spectral gap. 

When we add a next-nearest neighbor hopping t' to 
mimic imperfect nesting (as we expect to find generi- 
cally for realistic band structures, in particular in the 
face-centered cubic A3C60 materials), we find that the 
superconducting phase is only weakly affected but anti- 
ferromagnetism is strongly frustrated. As a result the su- 
perconducting region expands at expense of the magnetic 
insulator, and the superconducting dome may emerge in 
full [panel (b) of Fig. [5]. We propose that the grad- 
ual drop of superconducting order parameter past the 
dome maximum now nat urally explains the dec l ine of 
Tf. of expanded fullerides (iDahlke and Rosseinskvl l2002t 
iDurand oil . l2003t ICanin et aLl . l2008D . 

Finally, past the first-order Mott transition, we find 
that the antiferromagnetic insulator is formed mainly 
by spin-1/2 local configurations, which is in agreement 
with experiments in NH3K3C60 ( Iwasa and Takenobul . 
|2003( ). We also predict that ambient pressure A15 
CsaCgo, yet to be characterized, should similarly be a 
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FIG. 3 Superconducting solution (SC, left) and antiferromag- 
netic solution (AFM, right) anomalous self-energies (top) and 
order parameters (bottom) as a function of U/W. The top 
panels also show (green diamonds) the spectral gaps obtained 
multiplying A by the quasiparticle weight Z of each solution. 
The gaps are in units of the bandwidth W (the order param- 
eters are by definition dimensionless) . Dashed vertical lines 
mark the first-order phase transition between the two solu- 
tions. 



spin- 1/2 antiferromagnetic insulator. Besides spin rota- 
tional symmetry, this kind of state also breaks orbital 
rotational symmetry, signaling that spin ordering must 
be generally accompanied by orbital ordering. In am- 
moniated fuUerides that again is consistent with experi- 
ment (Iwasa and Takenobu, 2003;). We conclude that, in 
spite of strong simplifying assumptions, our model seems 
able to reproduce very important features of the phase 
diagram of expanded fuUerides. In the following we shall 
discuss in more detail the strongly correlated supercon- 
ducting phase near the Mott transition, and also propose 
experiments that might distinguish it from a standard 
BCS state. 



IV. UNDERSTANDING STRONGLY CORRELATED 
SUPERCONDUCTIVITY FROM DMFT AND THE 
IMPURITY MODEL 

The re-emergence of phonon-driven superconductivity 
close to the Mott transition - Strongly Correlated Su- 
percon ductivity (SCS) - was discussed in ICapone et all 
( 20021 ) in terms of Fermi-liquid theory. A key point of 
that phenomenon is the renormalization of the effective 
bandwidth and thus of the effective mass, both con- 
trolled (in a Bethe lattice) by the quasiparticle weight Z, 
Eq. (Uni). Z{U) decreases as a function of U /W and van- 
ishes at the continuous metal insulator transition point 
U — Uc, where the effective mass m*/m = 1/Z{U) di- 
verges, and the effective quasiparticle bandwidth W■^, — 
ZW vanishes. An estimate of the interaction between 
charged quasiparticles requires the evolution of fluctua- 
tions that take place in charge space. Because charge 



fluctuations are gradually frozen away near the Mott 
transition, the effective repulsion between quasiparticles 
is also renormalized down to some smaller value [/* < 
U. In particular, the Fermi-liquid description provided, 
e.g., by the Gutzwiller variational approach (Gutzwille^ . 
Il963( ) and supported by the DMFT behavior of the av- 
erage charge fluctuations {{n — 3)^), suggests that /7* ~ 
UTuZiU)^ - UZ{U), where Ff/ includ e s all the so- 
called vertex corrections ( Abrikosov et all . Il965l ). This 



implies that the vertex function Tu diverges close to the 
Mott tra nsition, but does no t compensate the vanish- 
ing of Z ( Capone et a/.l . 12^0021 ) . The pairwise Jahn- Teller 
and exchange based attraction J between quasiparticles, 
even if small, here is restricted to spin and orbital space, 
and has nothing to do with charge fluctuations. In other 
words, Hund's rule exchange and the JT coupling only 
influence the internal splitting of each molecular mul- 
tiplet without affecting its center of gravity. As a re- 
sult, this attraction should remain unrenormalized J* = 
JTj Z{U)^ ~ J close to the metal insulator transition, 
thus implying a strongly divergent vertex correction Tj 
that cancel the vanishing Z . Therefore the electron pair 
scattering amplitude in the Cooper channel should 
renormalize sis A = U - fj = Z{U) U ~ f J. 

When U jW is small, Z < 1, the main effect of U is to 
suppress superconductivity, as was noted earlier. How- 
ever, if \] is close to the critical metal-insulator value \Jc^ 
then Z ~ (C/c — U)/Uc <C 1 and the scattering amplitude 
turns negative in spite of a large U . This is qualita- 
tively the reason for the SCS re-entrance of supercon- 
ductivity (though in this region of course the actual pair 
scatte ring amplitude migh t deviate from this simple for- 
mula (jCapone et aZ.I . |2002| )). 



In addition, the Fermi-liquid argument suggests an ex- 
planation for the large value of superconducting order 
parameter, implying a large Tc, in the SCS regime, see 
Fig. [21 compared to the C/ = BCS values. In fact, 
when A^ ~ ZW, and Z is dropping sharply, the quasi- 
particle attraction A^^ will at some point U = U* equal 
the coherent quasiparticle bandwidth ZW . That very 
uncommon situation, of metallic quasiparticles with an 
pair attraction equal to their energy bandwidth, is known 
to yield maximum superconductivity for a given attrac- 
tio n. As shown by st udies of purely attractive mod- 
els (jMicnas et aZ.I . ll99Cll ). the maximum superconducting 
temperature ksTc attainable in that case is about 7% of 
the pair attraction energy itself. In our model of trivalent 
fuUerides, this estimate yields ksTc ^ 0.07A* ~ 0.2|J|, 
which has the correct magnitude of roughly 40 K for J ^ 
20 meV ~ a value in turn fully consistent with the ob- 
served spin gap 0.1 eV ~ 5J. While this coincidence of 
numbers is probably fortuitous, it does indicate that or- 
ders of magnitude implied by our model with realistic pa- 
rameters are quite consistent with experimental facts. At 
face value, it also suggests that 0.07 \ J\ — 0.07 (Jjt — Jh) 
could be the maximum attainable fc^Tc in fuUerides. We 
conclude that strong correlations play a crucial role in 
bringing the superconducting gap magnitude to the right 
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FIG. 4 Quasiparticle superconducting energy gap in units 
of the bandwidth W ~ 0.6 eV and on a larger scale than 
Fig. U computed through the anomalous self-energy, panel 
(a) in Fig. (2] multiplied by the quasiparticle residue Z(U). 
Notice that, above U /W ~ 1.1, the superconducting solution 
is metastable since the antiferromagnetic one has lower en- 
ergy, see Fig. (2] We note that the maximum gap Z Ac — 
0.045 W ~ 27 meV. 



range of values as compared with the experimental ones, 
see Fig. S) Such values and large critical temperatures 
would never be attained within conventional BCS theory 
using a value of A ~ 0.16 — 0.2, including as it should 
the large cancellation of JT by exchange. They could in 
point of fact be attained if the cancellation due to ex- 
change were (incorrectly) neglected; but then a lattice 
expansion should always lead to a decrease of Tc, con- 
trary to experiment. 

To appreciate further the effect of exchange- JT can- 
cellation, it is instruc tive to consider, as was done for 
a simplified model bv ICapone et al\ ( 2004I ). the behav- 
ior with U/W of the superconducting self-energy A(0) 
(proportional to the T = gap, and roughly speaking 
to Tc) starting with pure JT and without exchange, and 
then proceeding to turn on exchange and gradual can- 
cellation, see Fig. [5l For the bare JT, A ~ 1 (a strong 
coupling value) the superconducting self-energy decreases 
monotonically with increasing J7, in agreement with the 
Migdal-Eliashberg prediction of an increasing Coulomb 
pseudo-potential. Above a critical value, the system 
turns directly, via a second-order or weakly first-order 
phase transition, to a Mott insulating phase. This result 
is fully consistent with previous calculations by Han si 
et al. ( Han et al where the same type of Hub- 

bard model was studied within DMFT at finite temper- 
ature. Treating explicitly the electron-phonon coupling 
(including the full phonon dynamics) with A = 0.6 and 
neglecting exchange they obtained a superconductor with 
monotonically decreasing gap. 

Through a progressive reduction of A (mimicking JT 
cancellation by exchange) we find that a non-monotonic 
superconducting behavior makes its appearance as a 
function of U . Initially there is still a single super- 
conducting phase for all U/W values; but two differ- 
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FIG. 5 Anomalous self-energy A(0) (related to the supercon- 
ducting gap by a factor Z~^) for the superconducting solution 
and different values of the coupling parameter J = 0.02, 0.05, 
0.1 and 0.2, which correspond to A = 0.09, 0.21, 0.42, 0.85, 
respectively. Note that for large J, corresponding to a case 
where the Jahn Teller coupling is not canceled by Hund's 
rule exchange, superconductivity is strongest at (7 = 0. For 
increasing cancellation (decreasing J), two separate super- 
conducting pockets emerge, a BCS pocket near (7 = 0, and a 
SCS pocket near the Mott insulator. When the cancellation is 
nearly complete, SCS is many orders of magnitudes stronger 
than BCS. This is the situation we propose in our model of 
fu Uerides. Similar phy sics was described for a simpler model 
in lCapone et al\ \2QQ4 \. 



ent regions near zero and near Uc begin to materialize. 
(Note that Uc simultaneously shifts to higher C/ as A 
decreases). When the cancellation is so strong that A 
is still positive but small, the two superconducting re- 
gions break apart to form two separate pockets, leav- 
ing a normal metal phase in between. In the leftmost 
pocket near U/W = the anomalous self-energy has a 
BCS-like exponential dependence on A and indeed su- 
perconductivity in this corner is BCS. Superconductivity 
in the rightmost pocket near the metal insulator transi- 
tion behaves quite differently. Here the A dependence of 
superconductivity is much weaker, and the superconduc- 
tive gap much stronger, than in the BCS pocket. Super- 
conducti vity in this pocket c an in fact be characterized 
as SCS ( Capone et al\ . |2002| ). due to narrow quasipar- 
ticle pairing as described above. A similar behavior to 
Fig.O with two separate BCS and SCS regimes emerging 
from a single initial one when the effective pairing attrac- 
tion is progressively weakened by exchange, was derived 
and illustra ted in a simpler twofold degenerate model in 
ICapone et al. (2004). 

The SCS superconducting pocket near the Mott tran- 
sition is expected to differ from the BCS pocket even in 
its normal state properties. The normal state underly- 
ing a BCS superconductor is Fermi-liquid-like. On the 
other hand, previous analysis suggest that the Fermi- 
liquid picture is likely to break down in our model when 
the Mott transition is approached. The key reason for 
the breakdown of the Fermi- liquid is precisely that, when 
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Z 0, the attraction between quasiparticles must even- 
tually reach and exceed in magnitude the quasiparticle 
bandwidth ZW, a situation difficult to sustain.'^ Pos- 
sible deviations from th e Fermi-liquid paradigm were in 
fact overlooked in Rcf. ( Capone et all . [2002i) as they are 
related to the very-low energy behavior of the normal 
phase close to the Mott transition, not explored in that 
work. Later, the non-Fermi-liquid behavior was discov- 
ered in the two-orbital mod el where the physics is very 
similar ( Capone et a/.l . [2003 l. 



A. Anderson Impurity With a Rigid Bath 



The DMFT calculations described above involved two 
steps, one solving the Anderson impurity model (AIM), 
the other making that selfconsist ent with the ba t h. Fo l- 
lowing a reasoning proposed by iFabrizio et al. 1 (Hool), 
one may start off with the first step alone, namely analyz- 
ing the bare AIM without imposing any self-consistency 
constraint. The conduction bath can be assumed to 
have a flat density of states, and the bath-impurity hy- 
bridization to be structureless, a situation which avoids 
numerical uncertainties and yields accurate low-energy 
properties. This kind of analy s is app lied to the AIM 
© shows (|De Leo and Fabriziol . 120051 ) that two differ- 



ent impurity phases are stabilized according to the ra- 
tio betwe en the a. t tracti on J, and the Kondo tempera- 
ture Tk ( Hewsonl . Il997() . Below this temperature and 
when J = 0, the spin of an impurity coupled to a 
Ferr ni sea is screen ed out and absorbed in the conduction 
sea (|Hewsonl . [T997l ). In the lattice context within DMFT, 
the Kondo scale measures metallic coherence and corre- 
sponds to the renormalized quasiparticle bandwidth ZW . 
For finite J 7^ but smaller than Tk, Kondo screen- 
ing remains, thus still implying a Fermi-liquid behavior 
in DMFT. In fuUcrides, the impurity represents the Cgp" 
ion, carrying three orbitals and three spins. In the Kondo 
phase each of the three spins is separately screened by the 
bath and thus incorporated in the Fermi sea. Conversely, 
when J > Tk the Kondo screening is lost, and that was 
shown to imply a non Fermi-liquid phase characterized 
by a pseudogap in the single-part icle spectrum and by 
several other singular properties (i De Leo and Fabriziol . 
[1)05). A very qualitative description of this phase is 
that, unlike the Kondo phase, two spins out of three pair 
off antiferromagnetically at any given time, leaving out 
a single spin 1/2 available for Kondo screening. How- 
ever, since orbital degeneracy is unbroken, this residual 



^ Note that the high-energy Hubbard bands are unaffected by the 
small att raction J; supercond uctivity is just a matter of quasi- 
particles l lCapone et aLl . r2002[) . Therefore, a quasiparticle attrac- 
tion exceeding their bandwidth cannot correspond to an instabil- 
ity towards a Mott insulator, which must involve also the Hub- 
bard bands, but at most towards a breakdown of a quasiparticle- 
based Fermi liquid. 




FIG. 6 Impurity spectral function with a rigid bath charac- 
terized by a flat density of states (no DMFT self-consistency) 
across the critical point: the curves from top to bottom cor- 
respond to the evolution from the Kondo screened regime to- 
wards the pseudogap regime. The dotted line is the spec- 
tral function deep inside the non- Fermi-liquid phase J 3> Tk- 
The top inset shows the same curves plotted on a larger en- 
ergy range, where the pseudogap of the dotted line is vis- 
ible. The Fermi liquid behavior corresponds to the nor- 
mal state of the SCS superconducting phase for J < Jcp 
(analogous to the overdoped regime of cuprates) whereas the 
pseudogap behavior corresponds to the SCS normal state for 
J > Jcv (analogous t o the u nderdoped cuprates). Figure from 
iDe Leo and Fabrizid l|2005l ). 



spin = 1/2 also carries orbital momentum L = 1, 
which cor responds to an overscreened non-Fermi liquid 
situation Leo and Fabriziol . HqoI). In this regime, 
it was predicted that the impurity contributions to the 
specific-heat coefficient and to the pair susceptibility in 
the s-wave channel (|14p diverge as T~^/^ at low temper- 
ature T. In addition, the conduction electron scatter- 
ing rate has a non-analytic temperature-behavior T^/^. 
The local response functions to either a quadrupolar field, 
which splits the orbital degeneracy, or to a magnetic field 
which polarizes both spins and orbitals also diverge as 
y-i/s rpj^g phases, the Kondo screened phase and 
the pseudogap phase, are separated by a critical point 
at J = Jcp — Tk- It is endowed with a finite entropy 
1/2 In 3, and with a divergent superconducting suscepti- 
bility with an exponent 1/3. As shown in Fig. [HI the 
single-particle spectral function displays strong devia- 
tions from a normal metal in the pseudogap phase and 
at the critical point. 

Near this critical point it has been found that the low- 
energy dynamics around t he impurity is contro l led b y 
two separate energy scales ( De Leo and Fabriziol . l2005l) . 
T-|_ and T_ , whose behavior is very different as a function 
of increasing U/W. A higher energy scale T+ is set by 



the critical Jcp ^ Tk and represents the width of a broad 
incoherent resonance; it evolves smoothly and unevent- 
fully as the critical point is crossed (See Fig. jH]). A lower 
energy scale T- oc | J — Tk \ ^ measures instead the dis- 
tance from the critical point, and leads simultaneously to 



a narrow resonance in the Fermi-liquid region, and to an 
equally narrow spectr al density dip (the "ps eudogap" ) in 
the pseudogap phase ( De Leo and Fabrizid . [2005) . Since 
in this phase the impurity still carries a residual spin- 
1/2, there remains a finite value of the spectral function 
at the chemical potential, and the gap is not complete. 
The pseudogap widens if J is increased, and the cusp- 
like dip in the impurity spectral function smoothly turns 
into a cusp-like peak, the value at the chemical potential 
staying fixed and constant. This behavior is shown by 
the dotted curve in Fig. [5] corresponding to a very large 
pseudogap (see top inset), possessing a very tiny peak at 
the chemical potential. This indicates the existence of 
yet another energy scale besides T+ and T_ that sets the 
width of the cusp peak. 



B. Anderson Impurity with a Self-Consistent Bath in 
DMFT 

The rigid bath AIM behavior and its critical point 
briefly reviewed above provide a guide to the DMFT re- 
sults once the impurity-bath coupling is self-consistcntly 
determined. First of all, since Tr- coincides within DMFT 
with the renormalized ZW which in turn vanishes when 
the continuous Mott transition is approached, the impu- 
rity critical point is inevitably met before the metal insu- 
lator transition as U /W is increased, at some Ucp ^ Uc- 
This entails several important consequences: 

• The normal state may be a Fermi liquid only far be- 
low the continuous metal insulator transition, such 
as perhaps may be the case in the non-expanded 
fuUerides. Expanded compounds on the other hand 
are expected to have a non Fermi-liquid normal 
state, and eventually a pseudogap, possibly devel- 
oping before the first order transition to the anti- 
ferromagnetic insulator. 

• The SCS superconducting pocket near the Mott 
transition reflects the leading instability of the im- 
purity critical point. In other words SCS super- 
conductivity is the way in which the lattice model 
responds to impurity criticality and avoids it. 

• The lower energy scale vanishes right at the crit- 
ical point, r_ = 0. Here T+ ~ J thus remains 
as the only energy scale controlling the magnitude 
of the superconducting energy gap. Away from 
U = Ucp, T_ ^ 0, and the amplitude of the su- 
perconducting gap should decrease mon otonically 
with (T+ - T_)/T+ (jCapone et aLl . l2004D . since T_ 
cut-offs the local pairing instability. Therefore the 
gap should be maximum right at the impurity crit- 
ical point (the top of the dome). 

• Even though the normal phase is non-Fermi liquid, 
well defined Bogoliubov quasiparticlcs should exist 
inside the SCS superconducting pocket. 




u/w 



FIG. 7 Calculated weights of the zero- frequency delta- 
function contribution to the optical conductivity for the su- 
perconducting phase (superfluid stiffness, squares) and the 
normal phase (Drude weight, circles). The zero-frequency 
anomalous self-energy is also plotted. Note the higher values 
in the SCS superconductor, relative to the nonsuperconduct- 
ing metal phase. 



The last statement comes from the fact that, at the im- 
purity critical point, superconductivity provides a new 
screening channel which helps the system get rid of the 
finite residual entropy at the critica l point, thus eliminat- 
ing non-Fermi liqu id singularities ( De Leo and Fabrizid . 
120051 : ISchiro et all [2008) . 

The role of superconductivity as a novel screening 
channel close to the impurity model critical point should 
reflect, in the lattice model, into a gain of band energy 
(the tight binding "kinetic energy") at the onset of su- 
perconductivity. Owing to a sum rule connecting kinetic 
energy and zero-frequency optical con ductivity (some- 
times referred to as "Drude weight") (jScalapino et al\ . 
119931 ) the onset of SCS superconductivity close to the 
Mott transition should lead to a Drude weight increase. 
This prediction is well borne out by the full DMFT so- 
lution of our Hamiltonian. In Fig. [7] we plot the w = 
(d.c.) optical conductivity of our model superconductor 
(w here it coincides wit h the superfluid stiffness), defined 
by ([Toschi et a ll l2005f) 



kin 



X,,(q^O,r! = 0), 



(16) 



where Ekin is the kinetic energy and Xjj the static limit 
of the paramagnetic part of the electromagnetic kernel 



Xjj 



^ J deN{e)V{e) 



X [G(e,c^„)G*(e,c^„)+F(e,c^„)F(e,c^„)], (17) 



where V{e) 



(4^2 



^)/3 is 



the current vertex in 
the Bethe lattice, while G{e,ujn) and F{e,ujn) are the 
normal and anomalous lattice Green functions, respec- 
tively. In the same figure we also plot the zero-frequency 
d.c. conductivity (the Drude weight) of the underlying 
metastable solution where superconductivity is inhibited. 
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FIG. 8 Energetic balance underlying superconductivity. 
^Epot 



= E: 



pot — Epot is the diflterence between the poten- 
tial energies of the superconducting and the normal solution, 



while AEk 



is the same difference between 



the kinetic energies of the two solutions. 



this state meant to provide a cartoon of the real normal 
phase above Tc- The Drude weight is given by Eq. ([16]) 
and P7)l with F{e, w„) = 0. Upon entering the SCS dome 
from the low U/W side, the superconducting phase ini- 
tially loses kinetic energy over the normal state as in 
ordinary BCS theory. However, upon increasing U/W 
at and beyond the dome maximum, the loss reverts over 
to a gain, and indeed most of the SCS superconducting 
pocket is predicted to have a larger weight of the zero- 
frequency optical absorption than the non superconduct- 
ing state. The same behavior is displayed (Fig. [5]) by the 
energy balance between the two solutions. Only far be- 
low the Mott transition the superconductor is stabilized 
by a potential energy gain, as is the case in BCS theory. 
In the pseudogap regime, corresponding to the expanded 
fuUerides near the Mott transition, the stabilization is 
associated with a kinetic energy gain. 



A similar phenomenon is well known in the optica l 
cond uctivity of high-Tc copper oxides ( Carbone et all . 
I2006D . Our calculations show that an increase in the 



zero-frequency optical conductivity or a kinetic energy 
gain do not actually exclude an electron-phonon pairing 
mechanism, but rather demonstrates the key importance 
of strong electronic correlations. Thanks to pairing, the 
motion of carriers in the superconducting phase is facili- 
tated with respect to the pseudogap nonsuperconducting 
metal. In that anomalous metal - a non Fermi liquid- 
the interaction constraints jam the free propagation of 
quasiparticles, causing a kinetic energy cost, partly re- 
leased with the onset of superconductivity. It would be 
of extreme interest if the optical conductivity increase 
demonstrated in cuprates could be investigated in ful- 
lerides, both regular and expanded, since that would help 
discriminate between conventional BCS and SCS. 



V. DISCUSSION OF EXPERIMENTAL PROBES 

In the following we discuss whether and how the dom- 
inance of strong correlation we propose for the super- 
conducting fuUerides can be reconciled with the list of 
experimental facts given earlier, apparently supporting a 
conventional BCS behavior. The first quantity we dis- 
cuss is the specific heat jump at Tc- Within BCS theory, 
the specific heat jump AC'v at Tc is approximately 



ACv 



1.52 7* 



(18) 



where 7* is the specific heat coefficient of the metallic 
phase (Cy = 7*T)), proportional to the mass enhance- 
ment m* /m, in our case l/Z. The measured specific heat 
jump leads, through (fTS]) . to an estimate of 7* ~ 370. 
This is indeed a rather low and regular value which has 
been commonly advocated as evidence of weak correla- 
tions in fuUerenes. However, the standard argument only 
holds provided the normal phase is Fermi-liquid, which is 
not applicable close to a Mott transition, independently 
of any models. 

In Fig. [9] we plot the jump in Cy at Tc for increasing 
U/W, compared with the Fermi-liquid estimate 1/Z^ 
After a region where the calculated quantity closely fol- 
lows l/Z, ACv/Tc fiattens out and stays roughly con- 




UAV 

FIG. 9 Specific heat jump, in units of Tc, 70, as function of 
U/W. For comparison we also show the behavior of 1/Z{U), 
which should correspond to the same quantity if Fermi liquid 
theory were valid. 



* As mentioned in Sec. IIII.AI for the present three-orbital model 
it has not been practical to include in Eq. (|6)l a number of ex- 
cited states sufficient to make the truncation error in the Green's 
function negligible. For the data reported in Fig. |9]we have an 
average error of 3-4% which does not allow us to determine Tc 
with sufficient accuracy. Yet, the jump of the specific heat rela- 
tive to Tc turns out to be almost independent on the truncation 
error and on the details of the calculations. 
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FIG. 10 Normal-phase uniform magnetic susceptibility x 
normalized to the non- interacting value xo- For (7 = 0, 
X — Xo, the difference being extremely small since J ^ W . 
The plateau between U /W = 0.5 and 1 signals the effective 
crossover from a Fermi liquid with S= 3/2 per site, to a non 
Fermi liquid with 5=1/2 per site. 



slant around 470 up to the Mott transition, despite a 
diverging 1/Z (See Fig. [9]). This shows that the energy 
scale that controls the superconducting instability is con- 
stant near the Mott transition, consistent with the single- 
impurity prediction that this scale should be identified by 
T-|_, a quantity of order J. The bottom line conclusion 
here is that a normally sized specific-heat jump does not 
imply BCS superconductivity in fuUerides. 

Another physical quantity which seemingly pointed to- 
wards weak correlations in fullerides is the magnetic sus- 
ceptibility measured in the normal phase, apparently con- 
sistent with a weakly correlated a Fermi liquid and a 
Stoner enhancement of about a factor 2 to 3. This argu- 
ment again becomes inconclusive once the Fermi-liquid 
scenario is abandoned. In Fig. [TO] we plot as function of 
U /W the uniform magnetic susceptibility in the normal 
phase calculated by DMFT. After an initial Stoner en- 
hancement at small U /W , the susceptibility fiattens out 
and remains almost constant before a rapid growth which 
takes places extremely close to the Mott transition. In 
the plateau region the susceptibility enhancement of a 
factor between 2 and 3 with respect to U = Q, surpris- 
ingly close to the experimental enhancement, covering 
the whole superconducting domain. Physically, the ori- 
gin of this susceptibility plateau for increasing U is quite 
instructive. It corresponds to the gradual crossover of the 
maximum spin available at each site from S=3/2 in the 
Fermi liquid at small U < Ucp (Kondo screened AIM) to 
S=:l/2 near Mott and large U > Ucp, where Fermi liquid 
behavior is lost. In essence, at large U each molecule is ef- 
fectively in a (dynamically) JT distorted state, where two 
out o f three electrons are spin paired ( Auerbach et all . 
I1994D . We are thus led to conclude that the relatively 
weak observed enhancement of susceptibility does not 
correspond at all to a Stoner enhanced, weakly corre- 
lated Fermi liquid - in fact quasiparticles do not even 
exist in most of the plateau region. 

Finally, we wish to address signatures of the strongly 



correlated scenario which we expect should show up 
in important spectroscopies including the tunneling I- 
V characteristics and angle resolved photoemission spec- 
troscopy (ARPES). This is initially embarrassing on two 
accounts. First, ARPES spectroscopies are k-vector re- 
solved, whereas in DMFT we do not have access to any 
spatial structure. Second, tunneling spectroscopies are 
extremely well resolved near zero voltage, whereas our 
Lanczos method yields a much poorer spectral function 
resolution in this region. 

Let us addressing tunneling first. Although Fig. [H] 
refers to the impurity Cg^^ molecule, we believe that 
similar features would remain after full DMFT self- 
consistency in the normal phase - if we only had a bet- 
ter low-frequency numerical resolution than we presently 
have. Hence, we suggest that tunneling I-V spectra of 
expanded fullerides be measured and examined, in or- 
der to bring out the expected rich structure of the kind 
sketched in Fig. [HI 

Next, let us consider photoelectron spectroscopy. 
Agai n according to t he s ingle-impurity analy- 
sis ( De Leo and Fabriziol l2005t) . the imaginary part 
of the single-particle self-energy should be finite and of 
order T_|_ almost everywhere in the non-Fermi liquid nor- 
mal phase above Tc. This has the following implications 
for ARPES: 

• the fuUeride photoemission spectrum should show 
tiu bands dispersing in the Brillouin zone with 
nonzero bandwidth, governed by the energy scale 
T-f. The value of T4. decreases with increasing 
U/W (increasing expansion), from W at U = 
to (larger than) J at the Mott transition; 

• there should be a spectral peak broadening of 
the same order of magnitude T-(_ as the appar- 
ent k-resolved band dispersion. In particular the 
broadening should remain constant approaching 
the Fermi surface - unlike a Fermi liquid phase 
where quasi-particle peaks become narrower and 
narrower. 

The momentum-independence of the DMFT self- 
energy implies that, in this approximation, the k- 
modulation of the electronic dispersion is assumed to re- 
main unaffected by interactions. Within this assumption 
and without requiring too high frequency-accuracy, we 
can compute a toy k-resolved spectral function accord- 
ing to 



Aik,uj) 



-Im 



1 



cj — ek — S_dm_ft(w) ' 



(19) 



where ek is the non interacting dispersion and 
^dmft{^) is the DMFT self-energy calculated with a 
finite number of baths. The effect of the local self-energy 
will be to change the effective bandwidth and to give 
rise to finite lifetime effects, even if the k-modulation of 
the dispersion is unrenormalized. For our Bethe lattice. 
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FIG. 11 Simulated photoemission spectra elaborated form 
the DMFT result in the normal phase. The left panel refers to 
U/W = 1.1 (corresponding to unexpanded fuUerides), while 
the right panel is for U/W = 1.3 (corresponding to expanded 
fuUerides) . 



there is no straightforward definition of momentum, and 
we remedy that by computing A(e, uj), which corresponds 
to Eq. din]) with £k ^ £• 

In Fig.[TT]we show theoretical ARPES results for some 
choices of e obtained by using the DMFT self-energy for 
temperatures above Tc, both for a value of U/W — 1.1 
which lies close to the maximum of the superconducting 
dome, but still on the less correlated side, correspond- 
ing to unexpanded (or moderately expanded) fuUerides, 
and for a value which lies in the downward branch of the 
dome (U/W = 1.3), corresponding to a very expanded 
fuUeride. We note in both cases the existence of an inco- 
herent low-energy feature dispersing with a reduced but 
nonzero electron bandwidth of 0.1 W. In the expanded 
case the pseudogap feature is clearly present. 

Re cent photoemission spectra of K3C60 ( Goldonil . 
120071 ) indicate an overall dispersion bandwidth of about 
160 meV, about a quarter of the bare calculated band- 
width in the local-density approximation. The exper- 
imental spectral peak does not show the usual Fermi- 
liquid-like narrowing on approaching the Fermi level, a 
fact which is in agreement with our expectation for a 
non Fermi liquid (although a non expanded fulleride like 
K3C60 probably hes at the beginning of the SCS dome, 
where deviations from Fermi liquid are not massive). 
Experimentally the peak does not appear to cross the 
Fermi level, and the intensity instead drops, suggestive 
of a pseudogap. Unfortunately the spectrum shows very 
strong vibronic effects, reflecting the retarded strong elec- 
tron Jahn Teller coupling. This aspect is not covered by 
our unretarded approximation, but it heavily affects the 
line shape and hampers the extraction of purely elec- 
tronic features. Treatment of the vibronic effects, and a 
quantitative description of dispersion, will require aban- 
doning in the future our approximation of infinitely fast 
phonon dynamics, as well as a possible extension to clus- 
ter extensions of DMFT which allo w for different reno r- 
malizations of different momenta (|Hettler et a/.l . [2OO0I ) . 



VI. CONCLUSIONS 



Summarizing, we addressed the apparently contradic- 
tory properties of expanded trivalent fuUerides supercon- 
ductors and insulators - and to some extent of the whole 
family of fuUerides - and presented a theoretical scenario 
emphasizing the role of strong electron correlations. That 
is especially designed and appropriate for the more ex- 
panded members of the family, such as (NH3)3;NaK2C6o, 
LiaCeo, Cs3_a;Ka;C6o and Cs3_^Rb3;C6o, and the recently 
discovered A15 Cs^Ceo that are near or past the Mott 
transition. 

Our model explains the dome-shaped increase and sub- 
sequent decrease of Tc upon expansion of the lattice spac- 
ing in fuUerides; the coexistence of metallic behavior and 
of Mott insulator features such as the large NMR spin 
gap in all fuUerides, and the S ~ 1/2 spin in the in- 
sulating state (identified as a Mott- Jahn- Teller insula- 
tor). It explains why the s-wave Tc can be as high as 
40 K even though the Coulomb interaction strength is 
prohibitive, and why Tc does not automatically decrease 
upon increase of U /W. It also accounts for more stan- 
dard observations, such as regular specific heat jumps 
and moderately high spin susceptibilities, facts that were 
so far construed as evidence for conventional BCS super- 
conductivity. 

Besides those listed in the previous Section, one can 
anticipate a number of additional experiments that could 
provide "smoking gun" evidence for strongly correlated 
superconductivity in fuUerides. The tunneling TV char- 
acteristics observable, e.g., by a scanning tunneling spec- 
troscopy tip should, in an expanded fulleride, develop 
the low energy features typical of the Kondo impurity 
spectral function. The isotope effect upon carbon sub- 
stitution should also behave very unconventionally, and 
eventually get smaller as the superconducting dome is 
passed and the Mott transition is approached upon ex- 
pansion. In this regime, as the quasiparticle bandwidth 
ZW gradually falls below the typical energy Ticu of an 
increasing fraction of the eight Hg Jahn Teller modes, 
the associated retardation effect should in fact disap- 
pear. The expanded fuUerides and related materials, 
clearly not enough investigated so far, deserve in our 
view the strongest experimental attention. They com- 
bine elements that make them members at large of the 
high temperature superconductor family. They combine 
neighborhood of the Mott transition and predominance of 
strong electron correlations, with conventional elements 
such as electron-phonon s-wave pairing, that are typical 
of BCS systems. Our study identifies a pseudogap and 
other features in the IV tunneling spectrum, an increase 
of zero-frequency optical weight in the optical response 
of the superconducting phase, and the emergence of two 
separate energy scales in ARPES as the most urgent ex- 
perimental undertakings that could confirm of falsify our 
claims. 
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